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We investigate how Kondo phenomenon occurs in the Anderson model dynamically coupled with lo- 
cal Jahn-Teller phonons. It is found that the total angular moment composed of electron pseudo-spin and 
phonon angular moments is screened by conduction electrons. Namely, phonon degrees of freedom es- 
sentially contribute to the formation of singlet ground state. A characteristic temperature of the Kondo 
effect due to dynamical Jahn-Teller phonons is explained by an effective s-d Hamiltonian with anisotropic 
exchange interaction obtained from the Jahn-Teller-Anderson model in a non-adiabatic region. 
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Recent discovery of exotic heavy-fermion phenomenon 
in filled skutterudite compound SmC^SbiV' 2 has triggered 
rapid increase of renewed attention to quasi-Kondo phenom- 
ena with phononic origin. In particular, it has been clarified 
that the specific heat coefficient is almost independent of an 
applied magnetic field. 1 Namely, the heavy-fermion behav- 
ior in SmOs4Sbi2 is magnetically robust. It has been pointed 
out that this peculiar phenomenon can be understood by non- 
magnetic Kondo effect originating from phonons. 3 

In general, the Kondo-like phenomenon occurs in a con- 
duction electron system hybridized with a localized entity 
with internal degrees of freedom. When electrons are coupled 
with Einstein phonons, a double- well potential is formed in an 
adiabatic approximation, leading naturally to a two-level sys- 
tem, in which Kondo has first considered a possibility of non- 
magnetic Kondo behavior. 4 The two-level Kondo system ex- 
hibits the same behavior as the magnetic Kondo effect. 5 Fur- 
ther intensive efforts to include the low-lying levels of local 
phonon have been made to clarify the quasi-Kondo behavior 
in electron-phonon systems. 3 ' 6 

In the filled skutterudite structure, rare-earth ion is sur- 
rounded by the cage composed of twelve pnictogens. Then, 
the rare-earth ion easily moves around potential minima in 
off-center positions inside the pnictogen cage. This is called 
the rattling, which is considered to be a key ingredient with 
significant influence on electronic properties of filled skut- 
terudites. As a natural extension of the two-level Kondo prob- 
lem, Hattori et al. have analyzed the four- and six-level Kondo 
systems to consider the effect of rattling in filled skutteru- 
dites. 7 Their results seem to be consistent with the magneti- 
cally robust heavy-fermion behavior observed in SmOs 4 Sb 12 . 

Concerning the problem of rattling, the measurement of 
elastic constant of PrOs 4 Sb 12 has suggested that the off- 
center motion of Pr atom inside the Sb cage has degener- 
ate E g symmetry. 8 When we attempt to include the effect of 
off-center motion into a model Hamiltonian, a simple way 
is to regard the rattling as relative vibration of surround- 
ing atoms, leading to an electron-phonon interaction term. 
In a phenomenological level, the present author has consid- 
ered a linear coupling between degenerate /-electron orbitals 
with e u symmetry and dynamical Jahn-Teller phonons with 
E g symmetry in a multiorbital Anderson model constructed 
based on a j-j coupling scheme. 9-12 Numerical analysis of 
this model has revealed that quasi-Kondo behavior occurs due 



to the release of an entropy log 2 of the vibronic ground state, 
originating from the clockwise and anti-clockwise rotational 
modes of dynamical Jahn-Teller phonons. 9 This scenario is 
believed to have a potential to explain the Kondo effect with 
non-magnetic origin in filled skutterudites, 11 ' 12 but due to the 
complexity of the model including three orbitals, the mech- 
anism of quasi-Kondo phenomenon due to dynamical Jahn- 
Teller phonons has not been completely understood. 

In this letter, the Anderson model dynamically coupled 
with local Jahn-Teller phonons is numerically analyzed. Here 
spin degree of freedom is suppressed for simplicity, but orbital 
degree of freedom plays a role of pseudo-spin. It is shown 
that total angular moment J, composed of electron pseudo- 
spin and phonon angular moments, is screened by conduction 
electrons, leading to the singlet ground state of J = 0. A 
characteristic temperature Tk defined by a peak in the specific 
heat is discussed on the basis of the s-d Hamiltonian derived 
from the Jahn-Teller-Anderson model in the non-adiabatic re- 
gion. To complete the discussion, we also consider Tk in the 
adiabatic region. 

Let us introduce the Anderson model coupled with local 
Jahn-Teller phonons. In order to focus on the effect of cou- 
pling between electrons and local phonons, we consider ac- 
tive orbital degree of freedom, but for simplicity, spin degree 
of freedom is suppressed. Then, the model is expressed as 13 

H = Y^ £ kc{ T c k r + ^2(Vc f kT d T + h.c.) + H loc , (1) 

kr kr 

where Ek is the dispersion of conduction electron, Ck T is an 
annihilation operator of conduction electron with momentum 
k and orbital r, d T is an annihilation operator of localized 
electron with orbital r on an impurity site, and V is the hy- 
bridization between conduction and localized electrons. We 
set V = 0.25 and the energy unit is a half of the conduction 
bandwidth, D, which is set as 1 eV throughout this paper. 
The local electron term H\ oc is given by 

Hioc = U n a n b + pp + H cph , (2) 

where U denotes inter-orbital Coulomb interaction, n T — 
d\d T , p = n a + rib, and p is a chemical potential. In this 
paper, we consider the symmetric case with p = —{7/2. The 
electron-phonon coupling term iTph is given by 

H cph = g{Q 2 T x + Q 3 t z ) + (P| + P|)/2 

(3) 

+ lj 2 (QI + Q 2 3 )/2, 
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where g is the electron-phonon coupling constant, Q2 and Q3 
are normal coordinates for (x 2 — y 2 )- and (3z 2 — r 2 )-type 
Jahn-Teller phonons, respectively, Pi and P3 are correspond- 
ing canonical momenta, t x — S a db + d\d a , t z = d\d a — d\db, 
and ui is the frequency of local Jahn-Teller phonons. Note 
that the reduced mass is set as unity. When we introduce 
phonon operators a 2 and 03 through Q2 = (&2 + a 2 )/V^ 
and Q3 = (03 + c4)/ \f%j, i? C ph is rewritten as 



H cph = Lj^/a[{a 2 + a\)r x + (a 3 + a\)r z 
+ u)(a\a,2 + ala 3 + 1), 



(4) 



where a is the non-dimensional electron-phonon coupling 
constant, defined as a = g 2 /(2ui 3 ). 

Here we consider the conserved quantity of H^. 14 Using 
the commutation relation [r a , Tp] — 2ie aj a 7 r 7 with the Levi- 
Civita tensor e, we obtain [-ff C ph, J] — 0, where J denotes 
total angular moment, defined by J = L y + r y /2 with L y — 
— i(aga2— a\a 3 ) and t v — — i(S a db— dtd a ). It is convenientto 
introduce the electron and phonon bases to diagonalize r y and 
L y , respectively. Namely, we introduce b± — (a| ± ia\)/^/2 
and dj. = (d\ ± \d\)/ where a denotes a pseudospin and 
\ (J.) corresponds to +(— ) sign. Then, we obtain 



H loc = Un< t n l + fip + uj(b\b + + bib- + 1) 



+ uV2^[(b+ +bi)a + + {b- 

where n a = p — ri|+rij, cr+ - 

The total angular moment J is given by 

J = L z + a z /2, 



djcij,, and cr_ = u^up 



(5) 



dldr. 



d^df 



d\di 



(6) 
Cor- 



where L z = b + b+ — b_b- and a z 

responding to the above transformations, we also introduce 
c fcer = ( c ka ^ lc kb) I '"^2 for conduction electrons. Then, the 
whole Hamiltonian is still expressed as eq. (1), when we re- 
place orbital index r with pseudo-spin a. 

Let us briefly discuss the vibronic state of H\ oc at half- 
filling. The fcth eigen-state of H\ oc labeled by J is given by 



I* 



(fe,j) N 

loc / 



n=0 



[pi J ld\\J- 



V 



1 

2 ; '' 



(7) 



where \L;n) = \L + n,n) for L > and \n,n + \L\) for 
L < Owith |n+,n_)= {1/ ^n + \nA){a\ )"+ {al) n - |0) and 
the vacuum state |0) . The corresponding eigen-energy is given 
by E^ JS> + \l. Note that J takes half-odd-integer values as 
J = ±1/2, ±3/2, ±5/2, • • • and the coefficients satisfy the 
relations of p<Q 7) = and = pjjj. The double 

degeneracy in each eigen-state originates from the rotational 
mode of Jahn-Teller phonons. The ground state is specified by 
J = ±1/2 and k = 1. It is clearly observed that the phonon 
state with L z = J + 1/2 (J — 1/2) is coupled with the elec- 
tron state with pseudo-spin S z = —1/2 (+1/2), leading to 
the vibronic state specified by total angular momentum J. 

In order to investigate electronic and phononic properties 
of H at low temperatures, we usually discuss corresponding 
susceptibilities. The susceptibility of an arbitrary operator A 
is expressed by 

1 ^ e- E ~' T 



' " :,/T -\(m\A\n)\ 2 , (8) 
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Fig. 1. (Color online) (a) Coefficients in I^&q 1 '°') vs. a for U = 1 and 
ui = 0.5. (b) Coefficients in l^g 10 '} vs. n for a = 0, 1, and 5. 



where E n is the eigen-energy for the nth eigenstate \n) of H 
and Z is the partition function given by Z = ^ n e~ E "/ T . For 
the evaluation of susceptibilities, here we employ a numerical 
renormalization group (NRG) method, 15 in which momentum 
space is logarithmically discretized to include efficiently the 
conduction electrons near the Fermi energy and the conduc- 
tion electron states are characterized by "shell" labeled by N. 
The shell of N — denotes an impurity site described by the 
local Hamiltonian. The Hamiltonian is transformed into the 
recursion form as 



H 



N+l 



AH 



N +t N '^2{c < Ncr C N + lc r + C* N+la C Nc r), (9) 



where A is a parameter for logarithmic discretization, cn<t 
denotes the annihilation operator of conduction electron in the 
TV-shell, and t n indicates "hopping" of electron between iV- 
and (N + l)-shells, expressed by 

■A- X )(l- A-^- 1 ) 



(1 



2y/(l - A- 2JV - 1 )(1 

The initial term Hq is given by 



A -2JV-3) 



H = A- 1/2 [H, oc + Y,V(cLda 



dtcoa)}- 



(10) 



(11) 



We also evaluate entropy Si mp and specific heat C lmp of lo- 
calized electron. In the NRG calculations, a temperature T 
is defined as T = A _(A,_1 ^ 2 . In this paper, we set A=2.5 
and we keep M = 5000 low-energy states for each renormal- 
ization step. The phonon basis for each Jahn-Teller mode is 
truncated at a finite number Afp n , which is set as Af pn = 30. 

Before proceeding to the NRG results, it is useful to ex- 
amine the vibronic state of Hq. The fcth eigen-state of Hq is 
found to be expressed as 16 



I* 



(fe,J) v _' 



[^(cl t cl l + d\d\)\J;n) 



(12) 



E„ 



b i J l 

+ (4t4 _ rfT T c oi)l J ;"> 

where J takes an integer. Note that the ground state is charac- 
terized by J = and k = 1. 

In Fig. 1(a), we plot the coefficients af^, 6^°q, and cf^ vs. 
a in the ground state for U — 1 and u> = 0.5. At a = 0, the 
system is described by the Anderson model and we find the 
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singlet ground state in which total spin is zero. When we in- 
crease a, the component of cf l is rapidly increased. Namely, 
the phonon state with L z = +1 (—1) is coupled to the electron 
state with total pseudo-spin S^ ot = — 1 (+1), leading to the 
singlet state with J = 0. When we further increase a, multi- 
phonon components with n > 1 become significant, as shown 
in Fig. 1(b). The broad peaks in the coefficients are shifted to 
larger n with increasing a. In short, the singlet ground state 
with J = includes the significant component of phonon 
state with the moment L z = ±1, which is cancelled by the 
pseudo-spin moment of the electron state with 5* ot = =Fl. 
It is emphasized that phonon degrees of freedom essentially 
contribute to the formation of the singlet ground state, which 
is an important point of the Kondo effect in an electron system 
with dynamical Jahn-Teller impurity. 

Let us now discuss the NRG results of the Jahn-Teller- 
Anderson model. In Fig. 2(a), we show entropy and specific 
heat for a = 0.3, uj = 0.5, and U = 1. An entropy of 
log 2 originates from the double degeneracy in the local vi- 
bronic state with J — ±1/2, corresponding to clockwise 
and anti-clockwise rotational modes of dynamical Jahn-Teller 
phonons. 9 Here we define T* as a temperature at which the 
entropy becomes log 2. In Fig. 2(a), we obtain T* ~ 0.03. 
When we decrease temperature, the entropy of log 2 is even- 
tually released and we observe a peak in the specific heat, 
which defines a characteristic temperature Tk. In Fig. 2(a), 
we obtain Tk ~ 3 x 10~ 6 . The Kondo phenomenon due to dy- 
namical Jahn-Teller phonons is clearly observed for the case 
of Tk <C T*, which is satisfied in the non-adiabatic region. 

In Fig. 2(b), the results of several kinds of susceptibili- 
ties are shown. We find that Xo- x = X<j > Xa^ where 
o'x = + c_ and a y = — i(er + — <r_), due to the effect 
of Jahn-Teller phonons coupled with quadrupole operators r 2 
and r x , which are a x and a y , respectively, in the pseudo- 
spin notation. We remark that charge susceptibility immedi- 
ately vanishes at high temperatures, indicating that charge de- 
gree of freedom is irrelevant to the present Kondo-like phe- 
nomenon. On the other hand, pseudo-spin susceptibilities and 
X.j are suppressed around at T = Tk- Due to the discussion 
on the vibronic state of H , we deduce that the total moment 
J, composed of pseudo spin and phonon angular moments, is 
screened to form the singlet ground state of J = 0. In order 
to confirm that the present behavior is really the Kondo effect 
concerning J, we check that \J is scaled by a single param- 
eter Tk. In Fig. 2(c), we plot T\J as a function of T/Tk 
for a = 0.1 ~ 0.7. We actually find that T\j is expressed 
by a single function of T/Tk within the numerical precision, 
except for the high-temperature region in which the renormal- 
ization is not enough. Note that the pseudo-spin susceptibility 
does not satisfy such a scaling relation. 

In order to understand the behavior of Tk, it is useful to 
derive an effective model from H by using the perturbation 
expansion in terms of V. Note that the condition for the ex- 
pansion is given by uj ^> 2A, where A = TrpoV 2 is the width 
of the virtual bound state and po is the density of state of the 
conduction electron at the Fermi level. Namely, the following 
effective model is valid in the non-adiabatic region with large 
phonon energy. The adiabatic region will be discussed later. 

Here we distinguish the degenerate ground state of H\ oc as 



|R> = |d> 



(1,1/2) 
loc 



and |L) = K 



(l.-l/2)\ 



where "R" and "L" 
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Fig. 2. (Color online) (a) Entropy and specific heat for a = 0.3, ui = 0.5, 
and U = 1. (b) Susceptibilities for pseudo spin, charge, and total angular 
moment for the same parameters as in (a), (c) T\j as a function of the 
scaled temperature T/T K for a = 0.1 ~ 0.7. (d) T K /T^ vs. a for 
(7 = 1 and w = 0.5. Solid circles denote the numerical results, while a 
solid curve indicates the analytic result of the s-d model. 



denote the rotational directions and the ground-state energy is 
given 

b y s ioc 1/2) + M- After some algebraic calculations, we 



C kl C k'l)Sz 



obtain the effective s-d model as 

+ • 7 i(4a c fc'T 5 + + c U c k'iS-)}, 

where S z = (|R)(R| - |L)(L|)/2, S+ = |R)(L|, 
|L) (R|, and the exchange interactions are given by 



(13) 



2V 2 \p\ 



(l/2)l2 



Tf/n p (l,l/2) 

„=0 U/2- £ loc 



+ oj(2n + 1) 



and 



2V 2 [q\ 



(l/2)l2 



" 77 /9 7?( 1 ' 1 /2) 



S- = 



(14) 



(15) 



■w(2n + 2) 

Note here that J± > J z for a > 0. Namely, the significant 
effect of non-adiabatic phonons appears in the enhancement 
of the transverse component of the exchange interaction. This 
is consistent with \<j x = Xa y > X<j z for T > T K in Fig. 2(b), 
since Jahn-Teller phonons are coupled with a + and er_ . In 
other words, the directions of the rotational Jahn-Teller mode 
are easily converted for high-energy phonons. 

For the s-d model with anisotropic exchange interaction, 
Shiba has obtained the explicit expression for the binding 
energy E. 11 When we define the Kondo temperature Tk as 



T K 



-E, we obtain Tk as 



Tk = Z?exp 



-1 



2/W^ - Jl 



tan 



-i/v^T^f 



(16) 



In Fig. 2(d), we depict T K /T^ } vs. a, where t£ 0) is the 
Kondo temperature of the Anderson model with a = 0. Nu- 
merical results are shown by solid symbols. Note that Tk as 
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Fig. 3. (Color online) (a) Entropies for u = 0.02 ~ 0.5 with U = I 
and a = 0.5. Arrows denote the Kondo temperatures defined by the peaks 
in the specific heat, (b) T K /T^ vs. uj for U = 1 and a = 0.5. Solid 
symbol denotes the numerical results. Solid and broken curves indicate the 
analytic results for non-adiabatic and adiabatic cases, respectively. 



well as in the numerical results is defined as a tempera- 
ture which shows the peak in the specific heat. Analytic result 
is depicted by a solid curve, which indicates eq. (16) divided 
by = D e -i/(2JoPo) with J = W 2 /U. We find that the 
numerical results agree well with the analytic curve for the 
s-d model with the anisotropic exchange interaction. 

Thus far, we have concentrated on the non-adiabatic region, 
since we are interested in the Kondo effect due to dynamical 
Jahn-Teller phonons, but in order to complete the discussion, 
let us also discuss the adiabatic region. In Fig. 3(a), we show 
entropies for several values of oj with U = 1 and a = 0.5. 
When oj is decreased, T* is decreased and it gradually ap- 
proaches T^K In Fig. 3(b), T K /T^ is plotted as a function 
of oj for U = 1 and a = 0.5. For oj > 0.2, the numerical 
results are well reproduced by the analytic expression of the 
s-d model, since T* and are clearly separated. However, 
for smaller oj, the numerical results are significantly deviated 
from the solid curve of eq. (16). In the adiabatic region, as 
shown in Fig. 3(a), T* becomes comparable with T^' and 
the effective s-d model for non-adiabatic phonons loses its 
physical meaning. 

For adiabatic Jahn-Teller phonons, we intuitively under- 
stand the decrease of Tk with increasing a due to the effective 
enhancement of U by the static Jahn-Teller energy Sjt- 18-20 
In fact, the Kondo temperature in the adiabatic region is given 
by Tl d = L» e - 1 /(2Jo ff po) ? where j cff = W 2 /(U + 4E JT ) 

with E JT = aoj. In Fig. 3(b), we show T^/T^ ] by a bro- 
ken curve. It is considered that the numerical results should 
be gradually changed from the solid to the broken curve with 
decreasing oj. In order to observe such a gradual change more 
clearly, it is necessary to carry out the calculations for oj less 
than 0.01. However, in the present value of M = 5000, it is 
quite difficult to perform the NRG calculations for oj < 0.01 
with sufficient precision, since low-energy phonon excited 
states cannot be fully included. For the purpose to include 
the adiabatic phonons effectively, it is necessary to improve 
significantly the NRG calculations, not by simply increasing 
the value of M, For instance, it may be important to include 
phonon excitations as the polaron effect. We postpone such 
effort as one of future problems. 

Throughout this paper, we have ignored spin degree of free- 
dom. In spite of this simplification, we have found the inter- 
esting phenomenon that Jahn-Teller phonons contribute to the 
singlet formation. Note that this Kondo effect is related to or- 



bital degree of freedom. In actual materials, of course, there 
should exist active spin degree of freedom. In this sense, it 
is necessary to analyze the two-orbital Anderson model in- 
cluding Jahn-Teller phonons which are coupled with orbital 
(quadrupole) degree of freedom. In particular, it is an inter- 
esting problem to clarify how the Kondo effect due to dynam- 
ical Jahn-Teller phonons occurs in the spin-orbital system. It 
is another future problem. 

Since in general, a spin-orbit interaction is strong for / 
electrons, readers may think that it is meaningless to use the 
model with spin and/or orbital degrees of freedom for /- 
electron systems. However, by exploiting the j-j coupling 
scheme, it is possible to construct microscopic models for 
/-electron systems, 2123 in which spin degree of freedom is 
introduced to distinguish the degenerate states in the Kramers 
doublet, while orbital degree of freedom distinguishes the dif- 
ferent kinds of Kramers doublets. Thus, it is meaningful to 
consider the spin-orbital model for /-electron systems. 

Note, however, that relevant degree of freedom of / elec- 
tron is generally described by multipole. In the total angu- 
lar momentum J, t v is 2u octupole with xyz symmetry and 
L y denotes Jahn-Teller phonon angular moment, indicating 
that J does not include the magnetic moment. Thus, in prin- 
ciple, the Kondo phenomenon due to dynamical Jahn-Teller 
phonons is distinguished from the standard Kondo effect con- 
cerning the magnetic moment. However, the dominant mul- 
tipole degree of freedom is changed according to the crys- 
talline electric field (CEF) ground state. Furthermore, some 
multipoles are mixed depending on the crystal structure, CEF 
parameters, and /-electron number. The Kondo phenomenon 
due to dynamical Jahn-Teller phonons in actual /-electron 
materials such as Sm-based filled skutterudites will be dis- 
cussed in detail elsewhere. 12 

In summary, we have studied the Kondo phenomenon 
in the Anderson model coupled with dynamical Jahn-Teller 
phonons. It has been found that total angular momentum J 
composed of the electron pseudo-spin and phonon angular 
moments is screened by conduction electrons to form the sin- 
glet ground state of J = 0. We have shown that the char- 
acteristic temperature is well explained by the s-d model 
with anisotropic exchange interaction derived from the Jahn- 
Teller- Anderson model in the non-adiabatic region. 
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